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FOREWORD 

This  work  was  performed  in  whole  or  in  part  under 
the  combined  support  of  the  Department  of  the  Navy's 
Office  of  Naval  Research  and  Bureau  of  Weapons  [Contract 
Nonr  1841(40)],  the  Weapons  Guidance  Laboratory  of  the 
United  States  Air  Force's  Wright  Air  Development  Division 
[Contract  AF  33(6l6)-6046],  and  the  MITRE  Corporation 
under  Subcontract  No.  27. 
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ABSTRACT 

Correlation  techniques  can  be  applied  to 
infrared  detection  systems  to  improve  the  signal - 
to-noise  ratio  of  the  radiation  detector.  Two  such 
correlation  systems  are  compared  and  evaluated; 
(1)  a  single  detector  synchronous  rectification 
system,  and  (2)  a  two-detector  cross-correlation 
system.  The  two  systems  are  interpreted  in  the 
context  of  search  and  tracking  operations,  and  a 
discussion  of  the  transfer  of  target  information 
through  these  systems  is  presented.  An  experi¬ 
mental  study  verifying  the  theoretical  analysis  of 
the  two-detector  cross-correlation  system  has 
been  carried  out. 
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SYMBOLS 

Defined  by  Equation  (4,25) 

True  signal  component  contained  in  S 
True  signal  component  contained  in  s 
Reference  signal  waveform 
Frequency;  bandwidth  limits 
Target  signal  frequency 
Bandwidth  of  the  noise,  A  v(t) 

Equivalent  noise  bandwidth  (defined  by  Equation  2.20) 
rms  value  of  AN(t) 

Noise  component  contained  in  S 
Power  spectrum  of  Av(t) 

Maximum  value  of  P(f) 

Amplitude-normalized  power  spectrum  of  Av(t) 

Improvement  in  signal-to-noise  ratio  (defined 'by  Equation  4.11) 
Improvement  in  signal-to-noise  ratio  for  unit  bandwidth 
Unit-impulse  response  function 
Correlator  output  signal 
Defined  by  Equation  (4.8) 

Signal  from  radiation  detector  entering  correlator 

Integration  time 

Time 

Transfer  function  of  integrator  network 
A  constant  equal  to  1.781 

Root-mean-square  fractional  error  in  Q(T) 

Dirac  delta  function 

Amplitude  of  sinusoidal  target  signal 
Amplitude  of  sinusoidal  reference  signal 
Defined  on  pp.  8,  11,  and  12  with  various  subscripts 
rms  value  of  A  v(t) 

Noise  component  contained  in  s 
Time,  as  distinguished  from  t 
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SYMBOLS  (Concluded) 

4i(f)  Noise  power  spectrujn 

'I^t)  Autocorrelation  function 

(//(t)  Autocorrelation  function  of  A  v(t) 

'(t)  Autocorrelation  function  of  e'(t) 
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I.  INTRODUCTION 

In  the  past  decade  a  number  of  technical  reports  have  been 

written  on  the  uses  and  applications  of  the  correlation  function.  The 

foundation  for  this  development  was  laid  by  Wiener^  in  1930  with  the 

introduction  of  the  correlation  function  and  its  Fourier  transform,  the 

power  spectrum.  Applications  have  been  made  in  such  diverse  fields 

3  4 

as  astronomy  (Ref.  2,  pp.  120-126),  biophysics,  and  infrared  physics. 

One  of  the  most  promising  applications  of  the  correlation  function 

is  in  the  detection  of  signals  in  random  noise.  The  initial  papers  of  Lee 

5  6 

and  Wiesner  (1950)  and  Lee,  Cheatham,  and  Wiesner  (1950)  present 

a  discussion  of  the  autocorrelation  and  cross -correlation  methods  for 

the  detection  of  signals  in  noise.  A  somewhat  amplified  version  of  this 

7  8 

work  has  been  done  by  Fano  (1951)  and  Davenport  (1951).  A  mathe- 

9 

matical  analysis  of  the  Dicke  radiometer  together  with  another  system 
using  two  independent  receivers  has  been  carried  out  by  Goldstein^^ 
(1955).  Still  more  recently  Jacobson^  ^  ( 1 957)  has  examined  a  multiple 
receiver  system,  laying  primary  emphasis  upon  obtaining  directional 
information  from  the  received  signal.  In  addition  to  the  above,  the 
Massachusetts  Institute  of  Technology  (MIT)  Research  Laboratory  of 
Electronics  (RLE)  has  published  several  progress  reports  (e.g.,  Refs. 
12—14)  dealing  with  correlators  and  their  use  in  the  detection  of  signals 
below  the  noise  level. 

The  relationship  between  the  use  of  correlation  techniques  and 
optimum  filter  techniques  is  one  of  domain.  The  original  work  of 
North^^  (1943)  and  Van  Vleck  and  Middleton^  ^  (1946)  in  the  theory  of 

^Superscripts  refer  to  References 
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optimum  filters  is  reviewed  in  a  paper  by  Zadeh  and  Raggazzini 
(1952),  Also  appearing  in  this  paper  is  a  discussion  of  the  relationship 
between  correlation  techniques  and  optimum  filter  techniques.  Theo¬ 
retically,  no  claim  to  superiority  can  be  made  for  either  method, 
inasmuch  as  stationary  random  noise  contains  the  same  information 
whether  described  in  the  frequency  domain  (power  spectrum)  or  time 
domain  (autocorrelation  function).  However,  from  a  practical  point  of 
view,  the  operation  of  a  system  in  the  time  domain  may  be  more  feasible 
for  many  applications.  Comments  to  this  effect  may  be  found  in  Ref.  6. 

The  purpose  of  this  report  is  two-fold.  The  first  part  is  a  theo¬ 
retical  analysis  of  two  detection  s'ystems,  the  synchronous  rectification 
system  and  the  two-detector  cross-correlation  system.  The  function  of 
these  systems  is  to  improve  the  signal-to-noise  ratio  of  radiation  de¬ 
tectors.  The  analysis  followed  in  the  above  references  (viz.  5  —  11)  is 
generalized;  particular  attention  is  paid  to  the  effect  which  the  form  of 
the  power  spectrum  has  upon  the  improvement  in  the  signal-to-noise 
ratio.  Each  system  is  then  interpreted  in  the  context  of  the  current  in¬ 
frared  systems  being  used  for  search  and  tracking  operations.  The 
second  part  is  an  experimental  confirmation  of  the  two-detector  cross¬ 
correlation  system  by  means  of  a  correlator  developed  at  the  MIT  Naval 
Supersonic  Laboratory  (NSL). 
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II.  THE  SYNCHRONOUS  RECTIFICATION  SYSTEM 

A  block  diagram  of  this  system  is  shown  in  Figure  1.  The 
signal  and  noise  waveforms  included  in  the  diagram  indicate  the  suc¬ 
cessive  operations  which  improve  the  signal-to-noise  ratio  of  the  radi¬ 
ation  detector.  For  convenience,  the  target  signal  is  shown  as  a  sinusoid. 
Radiation  from  a  target  --  after  passing  through  the  chopper  --  falls  on  the 
detector  and  is  sent  to  the  preamplifier,  appearing  in  the  output  as  the  wave¬ 
form  indicated  in  Figure  1.  A  reference  signal  of  the  same  frequency  and 
phase  as  the  target  signal  is  produced  by  a  steady  light  source  and  de¬ 
tector  combination.  The  noise  in  the  reference -signal  waveform  is 
assumed  to  be  negligibly  small.  The  two  waveforms  are  subsequently 
multiplied  together  and  then  integrated;  the  result  being  that  the  target 
signal  will  increase  at  a  more  rapid  rate  with  time  than  the  root-mean- 
square  noise.  The  combined  operations  of  multiplication  and  integration 
(or  averaging)  are  called  c ross -correlation. 

The  system  indicated  in  Figure  1  may  have  several  alternative 
forms.  For  instance,  the  multiplier  can  be  replaced  by  a  phase-sensitive 
detector.  The  function  of  the  phase-sensitive  detector  is  simply  to  reverse 
the  sign  of  all  waveforms  during  the  negative  phase  of  the  reference  wave¬ 
form,  and  pass  all  waveforms  unchanged  during  the  positive  phase  of  the 
reference  waveform. 

Before  proceeding  with  the  analysis  a  few  definitions  must  be  made. 
The  output  of  the  radiation  detector  before  entering  the  correlator  will  be 
written  as 

s(t)  =  e  (t)  +  Av  (t)  ,  (2.1) 

wlicre  e(t)  and  Av  (t)  are  the  true  signal  and  noise  components,  respective¬ 
ly.  The  output  of  the  integrator  can  then  be  written  as 
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T 

0 


8(t)  e'  (t)dt 


(2.2) 


where  e'(t)  is  the  reference  signal  waveform  entering  the  correlator. 
Expanding  (2.2),  it  follows  that 


S  (t) 


e(t)e'(t)dt  + 


I, 


Av(t)  e'(t)dt 


(2.3) 


The  first  term  of  (2.3)  is  the  signal  component  in  the  output  of  the 
integrator. 


E(T) 


"I 


e(t)  e'(t)  dt 


(2.4) 


The  second  term  of  (2.3)  is  the  noise  component  in  the  output  of  the  in¬ 
tegrator. 


An  evaluation  of  the  second  term  of  (2.3)  requires  that  some 
of  the  properties  of  the  noise  Av(t)  be  known.  We  shall  assume  that 
the  cell  noise,  photon  noise,  and  circuit  noise  comprising  Av(t)  have 
the  following  general  properties: 


a)  Stationarity,  Av(t)*  =  Av(t±to)*  for  all  t© 


b)  <Av(t)2>.j.  =  <Av(t)2>.j, 


c)  Ergodicity,lim<Av(t)*  >_  =  Av(t)^ 

T— 00  ^ 

The  brackets  denote  a  time  average  over  the  period  T,  whereas  the 
bar  indicates  an  ensemble  (statistical)  average.  No  attempt  to  justi¬ 
fy  the  above  properties  will  be  made.  The  assumptions  are  not  unduly 
restrictive,  and  the  experimental  data  obtained  on  the  basis  of  these 
assumptions  substantially  confirm  their  correctness.  A  discussion 
of  the  above  properties  may  be  found  in  Reference  18  together  with 
appropriate  source  material. 

It  is  well  known  that  the  output  of  a  linear  network  can  be  ex¬ 
pressed  by  the  equation 
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AN(t) 


An{T)Ro  (t  —  t)  dr 


(2.5) 


where  An(T)  is  the  input  waveform  to  the  linear  network  and 
Rq  (t  —  t)  is  the  unit-impulse  response  function  of  the  network. 
This  equation  can  be  verified  by  a  simple  application  of  the 
Fourier  integral  theorem.  To  find  the  statistical  mean-square 
value  of  AN(t)  we  make  the  following  operations: 


AN(t)2 


An(T)Ro(t  —  t)  dr  \  An(T')Ro  (t  —  t')  dr' 


(2.6) 


p  00  p  00 

AN(t)^  =  \  \  An(T)An(T’)J^  (t  —  t)Ro  (t  —  t')  dr  dr*  ,  (2.7) 

*^-00  ^  -00 


where  t'  is  introduced  to  distinguish  the  two  independent  integra¬ 
tions.  Taking  the  statistical  average  of  (2,7)  and  making  use  of 
property  b)  above,  we  have 


AN(t)^ 


An(T)An(T')  Ro  (t  —  t)Ro  (t  —  t')  dr  dr' 


(2.8) 


The  stationarity  property  permits  us  to  equate  An(T)An(T')  with  the 
autr covariance  function  —  of  the  noise  An(T);  that  is, 


AN(t? 


'^{r  —  t')Ro  (t  —  t)Ro  (t  —  T')dT  dr' 


(2.9) 


-00  ^  -00 


The  ergodic  prope  'ty  implies  the  equality  of  the  autocovariance  and 
autocorrelation  functions.  Moreover,  since  the  power  spectrum  of 
a  random  process  is  defined  as  the  Fourier  transform  of  the  autocor¬ 
relation  function,  we  may  write 


-T'l 


$  (f)  cos  2lTf(T  —  t')  df 


(2.10) 


Equations  (2.9)  and  (2.10)  require  only  that  the  power  spectrum  of 
the  noise  and  the  unit-impulse  response  function  be  known  in  order 


TR  307 


CONFIDENTIAL 


CONFIDENTIAL 


to  obtain  the  mean-square  noise  in  the  output  of  the  network. 

The  autocovariance  function  of  a  stationary  random  process, 
e.g.  u(t),  is  defined  by  the  equation 


4'{t)  =  u(t)  u  (t  +  t)  (2.11) 

where  the  bar  denotes  an  ensemble  average.  Some  of  the  properties  of 
the  autocovariance  function  include  the  following: 

'I'(t)  =  -  t)  (2.12) 


and 


^(t)  <  ^(0)  =  u^ 


(2.13) 


where  u^  is  the  ensemble  mean-square  value  of  u(t).  The  autocor¬ 
relation  function  is  defined  by 

1 

'k(T)  =  lim  yrrr  \  u(t)u(t  +  T)dt  (2.14) 

T^OO.  J-T 


In  view  of  the  ergodic  property  already  assumed,  the  autocorrelation 
and  autocovariance  functions  are  equivalent. 

The  unit-impulse  response  function  of  the  integrator  network  of 


(2.3)  is 


Ro  (t-T) 


6  (t— T)dt 


(2.15) 


where  6(t  —  t)  is  the  Dirac  delta  function.  Equation  (2.15)  is  very 
easily  evaluated  to  give 
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fl  ,  0  <  T  <  T 
R®  ,T>TorT<0 


(^.16) 


Inserting  (Z.IO)  and  (2.16)  into  (2.9)  and  carrying  out  the  two  inte¬ 
grals  in  T  and  t'  gives 

=  1  r“  $(f)(l  -co8ZTrfT)df  .  (2.17) 

^  ^Loo  £2 

The  form  of  the  power  spectrum  $  (f)  depends  upon  the  power  spectrums 
(or  autocorrelation  functions)  of  the  signals  entering  the  multiplier  and 
the  transfer  function  of  the  multiplier. 

Since  the  two  functions  in  the  second  term  of  (2.3)  are  inde¬ 
pendent,  the  autocorrelation  function  of  the  product  Av(t)e'(t)  can  be 
written  as 


’i'(T)  =  »|j(t)  (t) 


(2.18) 


where  0  (t)  and  0'(t)  are  the  autocorrelation  functions  of  Av(t)  ande'(t), 
respectively.  Using  the  fundamental  theorem  that  the  autocorrelation 
function  and  the  power  spectrum  constitute  a  Fourier  transform  pair 
(see  Eq.'(2. 1 0)),the  power  spectrum  associated  with  ^(t)  can  be  written  as 


•  r.a 


P(f)  cos  ZirfT  df  Fil;'  (t)  cos  ZirfrdT 


(2.19) 


where  P(f)  is  the  power  spectrum  of  Av(t).  For  convenience  a  noise 
equivalent  bandwidth  will  be  defined  as  follows: 


P(f)  df 


(2.20) 


where 


P(f) 

"max 


In  addition  the  mean-square  value  of  Av(t)  will  be  written  as 

pOO 

=  \  P(f)df  =  P  Af”''*  . 

J  ^  max 


(2.21) 


(2.22) 


Inserting  (2.19)  into  (2.17)  and  using  the  results  of  (2.20)  and  (2.22), 
the  mean  square  noise  in  the  output  of  the  integrator  becomes 
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AN^=—G{T)  ,  (-i.23) 

where 


Making  the  definition  hP  =  AN^  ,  the  root-mean- square  noise  becomes 


The  general  equation  for  the  signal-to-noise  ratio  obtained 
by  synchronous  rectification  can  now  be  written.  Taking  the  ratio 
of  (2.4)  to  (2.25),  we  have 

T 

E  1  e(t)e'(t)dt  ^ 

This  equation  leads  to  the  following  conclusions.  The  signal-to- 
noise  ratio  obtained  by  the  synchronous  rectification  system  is 
dependent  upon  three  functions: 

1)  the  form  and  frequency  of  the  target  signal 

2)  the  form  of  the  reference  signal 

3)  the  power  spectrum  of  the  noise  entering 
the  correlator 

A  more  detailed  discussion  of  the  effect  of  these  three  functions  upon 
(2.26)  will  be  given  in  Section  IV. 

In  applying  the  analysis  to  a  system  having  real  components, 
deviations  from  the  ideal  transfer  functions  of  the  multiplier  and 
integrator  must  be  included  in  (2.Z6).  This  means  that  the  Fourier 
series  coefficients  of  e(t)  and  e'(t),  and  the  power  spectrum  p(f) 
will  take  on  a  somewhat  modified  form. 
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III.  THE  TWO-DETECTOR  CROSS-CORRELATION  SYSTEM 


A  block  diagram  of  this  system  is  shown  in  Figure  2.  The  sig¬ 
nal  and  noise  waveforms  included  in  the  diagram  again  indicate  the  suc¬ 
cessive  operations  which  improve  the  signal-to-noise  ratio  of  the  radi¬ 
ation  detector.  Radiation  from  a  target  passes  through  the  chopper  and 
falls  on  both  of  the  detectors.  The  signals  are  then  amplified  and  ap¬ 
pear  in  the  outputs  of  the  pre-amplifiers  as  the  waveforms  indicated  in 
Figure  2.  The  signals,  containing  both  tme  signal  and  noise  components, 
are  then  multiplied  together  and  integrated. 

The  analysis  of  this  system  follows  the  same  pattern  as  the  pre¬ 
ceding  one.  The  two  signals  entering  the  multiplier  will  be  written  as 


Si  (t)  =  e^  (t)  +  Avi  (t),  (3.1 ) 

and 

S2(t)  =  e2(t)  +  Av2(t),  (3.2) 

where  e.(t)  and  Av.(t)  represent  the  true  signal  and  noise  components, 
respectively.  The  output  of  the  integrator  can  then  be  written  as 


S(t) 


=  J  Si  (t)s2 


(t)dt 


(3.3) 


or 


S( 


t)=  ^  [e,  (t)e2(t)  +  e,  (t)Av2(t)  +  e2(t)Avi  (t)  +  Av,  (t)Av2(t)]dt  .  (3,4) 


The  first  term  of  (3.4)  is  the  signal  component  in  the  output  of  the 
integrator, 
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C  C  i'\i  r  i  D  L  N  f  !  A  L 


E(t) 


(t)e2  (t)dt 


(3.5) 


The  last  three  terms  of  (3.4)  are  the  noise  components  in  the  output 
of  the  integrator. 

The  noise  from  these  three  terms  can  be  written  as 

AN(t)=  AN,  (t)  +  AN^Ct)  +  AN3{t)  .  (3.6) 

where  AN,  (t).  AN,(t).  and  AN3(t)  are  defined  by  the  equations 


and 


an,  (t)=j  e2(t)Av,  (t)dt 
0 

X 

AN2(t)=  r  e,(t)Av2(f  ),ll 
'^0 

AN3(t)sJ  Av,  {t)Av2(t)dt 


(3.7) 


(3.8) 


(3.9) 


Since  the  above  three  random 
can  easily  be  shown  that 


noise  variables  are  independent,  it 


aTF  =  ANj?  +  +  ANj^ 

The  first  two  terms  in  (3  101  havp  i, 

m  nave  already  been  determined  (see 

(<i.3),  (2.23),  and  (2.24)).  The  result  is  that 


(3.10) 


where 


and 


N,^=^0i(T) 

N,2  H  AN,^  ,  V  2  =  ^ 


(3.11) 
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(3.12) 


Similarly, 


N|=j^8,(T)  , 


(3.13) 


where 


0*(T) 


(tIcos  Zirfr 


[!' 


P2  (f)coB  Zirfr  df 


( 1  —  cos  2wfT) 
- p - 


df 


poo 

\  Pi(f)df 
-00 


(3.14) 


The  autocorrelation  functions  ^i(t)  and  ^2(t)  are  those  corresponding 
to  Cj  (t)  and  02  (t);  Pj  (f )  and  P2  (f)  are  the  amplitude  normalized 
power  spectrums  corresponding  to  Avj  (t)  and  Av2(t). 

The  last  term  of  (3.10)  can  be  found  by  altering  the  form  of  (2.23) 
and  (2.24).  If  the  function  is  redefined  as  the  autocorrelation  func¬ 
tion  of  Av2(t),  then  can  be  written  as 

p  00 

l/' '(t)  =  \  P2  (f)  cos  2’irfT  df  ,  (3.15) 

-00 


where  P2(f)  is  the  power  spectrum  of  Av2(t).  It  therefore  follows  that 

N3"=^^03(T)  ,  (3.16) 


where 


flj (T) s 


C  -I  r  ^  Pl  (f)cos  ZnfT  df  r  pj  (f)cos  2 Tifr  df  coB2TrfTdTl  — - 

s  "  PC  -  00  -  00 _ *•'  "00 _  _ J _ 


—  cos  ZwfT) 


df 


roc  ^00 

p,  (f)df\  pj(f)df 


(3.17) 
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Using  (3.11),  (3.13),  and  (3.16)  the  root-mean-square  noise  in  the  out¬ 
put  of  the  integration  can  be  obtained: 


(3.18) 


Taking  the  ratio  of  (3.5)  to  (3.18),  there  follows  that 


^  e,  (t)e2(t)dt 


^FeVfl 


t 


(3.19) 


where 

0(T)=  03(T)  +  -^  02(T)  +;^0,(T)  .  (3.20) 

This  equation  should  be  compared  with  (2.26).  The  signal-to-noise 
ratio  obtained  by  cross -correlation  is  dependent  upon  the  form  and 
frequency  of  the  target  signals  and  the  power  spectrums  of  the  noise 
entering  the  correlator.  In  the  following  section  a  detailed  expres¬ 
sion  for  the  ratio  E/N  will  be  discussed.  As  in  the  previous  system, 
deviations  from  the  ideal  transfer  functions  of  the  multiplier  and  in¬ 
tegrator  must  be  included  in  (3.19),  if  the  analysis  is  to  be  inter¬ 
preted  in  terms  of  real  system  components. 
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IV  .  IMPROVEMENTS  IN  SIGNAL- TO-NOISE  RATIO 

The  results  of  Sections  II  and  III  provide  a  general  description 
of  the  improvement  in  the  signal-to-noise  ratio  which  can  be  obtained 
with  each  system.  A  more  detailed  description  requires  that  assump¬ 
tions  be  made  concerning  the  form  and  frequency  of  the  target  and 
reference  signals  as  well  as  the  form  of  the  power  spectrum  of  the  noise 
entering  the  correlator.  In  the  analysis  which  follows  the  synchronous 
rectification  system  will  be  treated  first. 

The  true  signal  component  of  (2.1)  we  shall  assume  to  be  a  sine 
wave  of  frequency  fo  and  amplitude  c  , 

e(t)  =  e  sinZirfot  .  (4.1) 

The  reference  signal  will  be  taken  as 

e'(t)  =  c '  sin  Zirfot  .  (4.2) 

This  equation  permits  the  determination  of  the  autocorrelation  function 
\Ij  '(t)  appearing  in  (2.24).  From  (2.14)  it  follows  that 

^ '(t)  =  cos  2TrfoT  .  (4.3) 

Furthermore,  the  integral  appearing  in  (2.2  6)  can  now  be  evaluated. 
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For  fo  T  >  10  (  <  1%  error), 

E  =  — ^  T  .  (4.5) 

The  remaining  function  needed  to  evaluate  (Z.Zb)  is  the  power 
spectrum  of  the  noise  entering  the  correlator.  Some  difficulty  is 
encountered  in  selecting  a  power  spectrum  which  is  sufficiently 
general  in  form  and  yet  reasonable  to  handle  mathematically. 

Perhaps  the  most  satisfactory  one  has  the  following  form: 


r  1  ,  for  f,  <  f  <  fj  ,  -  f^  <  f  <  -  f, 

^  0  ,  fo  r  —  fj  <  f  <  f ,  ,  f  >  f  2  ,  f  <  —  f  2 


(4.6) 


where  f|  and  f2  determine  the  bandwidth  Af  of  the  noise.  A  diagram 
of (4.6) is  given  in  Figure  3. 

On  the  basis  of  (4,3)  and  (4.6),  the  function  0(T)  of  {Z.Zb)  may 
be  evaluated.  The  result  is  that 


e(T)  = 


1111 
+  fi  )  f2  )  ZHo  -  f,  )  2  (^2  -  fo  ) 


(4.7) 


cos  2-trT(fo  —  £|  )  ,  cos  2TrT(f,  —  fn  )  .  cos  2TrT(fo  +  f,  ) 
2(fo-fi)  Z(l2  ^ 


cos  2 
ZJT, 


+  ttT  Si  2TiT(fo  +  fj  )  -  ttT  Si  2uT(fo  +  f,  ) 
ttT  Si  2TTT(fo  -  f,  )  +  ttT  Si  2TTT(f2-  f o  )  J- 


where  the  function  Si(x)  is  defined  as 

Si(3 


/^X 

-/  V  —  I  sinu  j 
)i(x)  =  \  — —  du 


(4.8) 
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The  above  function  converges  very  rapidly  to  u/Z  as  x  oo;  for  x  >  10 
Si{x)  s  tt/Z  (witliin  6  %).  An  exhaustive  treatment  of  the  function  can 
be  found  in  Ref.  19,  pp.  1-9. 

The  Equation  (Z.Z6)  can  now  be  written  explicitly: 


E 

N 


sin  4n  fp  T  \  T 

/Kje{T) 


(4.9) 


where  0(T)  is  given  by  (4.7).  The  signal-to-noise  ratio  obtained 
by  synchronous  rectification  depends  upon  a  complex  function  of  the 
chopping  frequency,  bandwidth,  and  integration  time. 

A  close  study  of  (4.9)  reveals  that  a  number  of  important  trends 
can  be  determined.  E/N  approaches  zero  as  T^  for  very  small  T, 
and  varies  as  the  square  root  of  the  integration  time  for  very  large  T. 
Between  the  limiting  values  of  T  oscillations  occur  in  a  somewhat 
complex  manner. 

The  effect  of  the  chopping  frequency  fo  and  bandwidth  Af  upon  the 
improvement  in  the  signal-to-noise  ratio  is  also  significant.  The  most 
satisfactory  way  to  determine  this  effect  is  to  discuss  the  form  of  the 
power  spectrum  of  the  multiplied  noise  (see  (Z.19)).  The  reason  for 
this  approach  is  that  the  integrator  network  following  the  multiplier 
operates  as  a  low-pass  filter.  The  transfer  function  of  an  integrator 
network  integrating  over  a  period  of  time  T  can  be  shown  to  have  the 
following  form: 

‘Wi(f)  =2^  2wfT  +  j(cos  ZufT—  ijj'  . 


In  terms  of  this  equation  (Z.17)  becomes 

N*  =  f'”§(f)  |Wj{f)|2  df  .  (4.10) 

-  00 

The  function  |  W.(f)|2  is  shown  in  graphical  form  in  Figure  4. 
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It  is  apparent  from  the  above  two  equations  (and  from  that 

the  mean-square  noise  in  the  output  of  the  integrator  is  much  smaller 
for  power  spectrums  having  components  predominantly  at  high  frequency. 
We  now  proceed  to  a  discussion  of  the  power  spectrum  of  the  multiplied 
nois  e . 

When  (ii.l9)  is  evaluated  three  cases  are  observed  in  the  result: 

(a)  fo  -  fi  <  fz  -  fo  <  fo  +  ;  (b)  fz  -  fo  <  fo  -  fi  <  fo  +  ;  and  (c)  fo  -  fj  < 

fo  +  fi  <  fz  —  fo  •  The  power  spectrums  corresponding  to  these  cases  are 
shown  in  Figure  5  a,  b,  and  c.  These  graphs  show  that  low  frequency 
components  appear  in  the  power  spectrum  of  the  multiplied  noise  that 
did  not  appear  at  all  in  the  power  spectrum  of  the  cell  noise  (see  (4.6)). 
Furthermore,  no  manipulation  of  the  bandwidth  or  signal  frequency  fo 
can  eliminate  these  low  frequencies.  The  most  that  can  be  done  is  to 
minimize  the  contribution  to  the  integral  in  (4.10)  due  to  the  low  fre¬ 
quencies,  and  push  as  much  of  the  power  spectrum  as  possible  into  the 
high  frequency  range.  To  accomplish  this  fo  should  be  set  equal  to  f 3  ; 
the  chopping  frequency  fg  should  be  as  high  as  possible,  since  (from 
Figure  5  a,  b)  the  high  frequency  components  in  the  power  spectrum  are 
dependent  upon  its  magnitude.  The  results  of  prescribing  these  condi¬ 
tions  is  that  the  value  of  the  integral  in  (4.10)  will  be  minimized. 

Having  discussed  qualitatively  the  effects  of  integration  time  and 
power  spectrum  upon  the  root  mean-square  correlated  noise,  (4.9)  may 
now  be  represented  graphically.  However,  let  us  first  define  a 
quantity  Q(T)  as  follows: 


Q(T)  - 


(E/N) 

(e/v) 


(4.11) 


Q(T)  represents  the  improvement  in  signal-to-noise  ratio  obtained  for 
a  given  integration  time  and  noise  bandwidth.  Written  in  terms  of 
Q(T),  (4.9)  becomes 


Q(T)  e 

42 


sin4TTfoT  \  T 

/  4wn 


(4.1Z) 


TR  307 


16 


CONFIDENTIAL 


CONFIDENTIAL 


Figure  6  shows  a  graph  of  (1.12)  for  Q  (T)  >  1  and  for  various  values 
of  the  parameter  A  f .  The  signal  frequency  fg  has  been  set  equal  to 
the  high  frequency  limit  £2  of  the  band  pass.  The  function  Q(T),  as 
shown  in  the  graph,  increases  as  the  square  root  of  the  time -bandwidth 
product.  It  should  be  pointed  out  that  the  dependence  of  Q(T)  upon  the 
square  root  of  the  bandwidth  Af  is  somev/hat  misleading,  since  the 
initial  signal -to-noise  ratio  is  inversely  proportional  to  the  square 
root  of  the  bandwidth. 

A  graph  of  an  additional  function  Q’'‘(T)  is  shown  in  Figure  7. 

This  function  represents  the  improvement  in  signal-to-noise  ratio  for 
noise  of  unit  bandwidth  (£2  —  f  1  =  1).  The  band  has  been  moved  over 
the  frequency  spectrum  by  assigning  values  to  the  center  frequency 
f^  ^  (fj  +  £2)12..  We  may  conclude  from  the  figure  that  for  improve¬ 
ments  of  magnitude  Q(T)  >  1,  the  system  operates  nearly  independently 
of  where  the  bandwidth  Af  is  located  in  the  frequency  spectrum.  In 
other  words,  the  advantage  obtained  with  high  chopping  frequencies  is 
negligibly  small. 

Proceeding  to  the  two-detector  cross -correlation  system,  we 
shall  make  the  same  assumptions  with  regard  to  the  target  signals  and 
the  power  spectrums  of  the  cell  noise  as  in  the  previous  system.  The 
true  signal  components  in  (3.1)  and  (3.Z)  will  be  written  as 

e|  (t)  =  e  1  sin  Zirfo  t 
^2(1)  =  Ee  sinZirfot 

The  autocorrelation  functions  appearing  in  0i(T)  and  02  (T)  (see 
(3. 12)  and  (3.14)  can  be  obtained  from  (2.14): 

(t)  =  cos  ZTTfgT  ,  (4.15) 

and 

^^2  (t)  =  ^-2—  cos  Zirfo  T  •  (4.16) 

Using  the  above  two  equations,  (4.13)  and  (4.14),  the  integral  in 
(3.19)  can  be  evaluated  to  give 
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j’^e,(t)e,(t)d.=4£i  .  (4.17) 

0 

For  fQ  T  >10  (<  one  per  cent  error), 

T 

J  ei(t)e2(t)dt  =  T  .  (4.18) 

The  power  spectrums  of  the  cell  noise,  Pj  (f)  and  PzCf),  are  assumed 
to  be  identical,  and  of  the  same  form  as  (4.6);  that  is. 


Pi  (f)  =  Pz  (f)  = 


Lo  , 


for  £i  <  f  <  iz  ,  —  fz  S  ~ 
for  —  f|  <  f  <  fj  ,  f  >  f2  ,  f  <  —  fj 
With  the  above  assumptions  the  equation  (3.19)  can  be  carried  out 
explicitly.  The  result  is 


(4.19) 


where 


E  _  tr  /  1  _  sin4iTf  0  T  \  £j_  T 


(4.20) 


0(T)  =  03  (T)  +  —  02(T)  + -J— 0i(T)  ,  (4.21) 

Vi^ 

and 

03  (T)  = - - - i-Z  -  2  1ny  -  In4f,f2  +  21n  j|2-±iL>  -  21n2TTT 

Mfz  -  fi)"  I 


+  cos  4iTf2  T  +  cos  4TTfi  T  +  4  sin  2TTf2  T  sin  2TTfj  T  +  Ci(4TTf2  T) 


+  Ci(4Trf,  T)  -  2Ci^Tr(f2  +  fj  )T  +  2Ci  2TT(f2  -  f,  )T  +  4T7f2  T  Si  4irf2  T 
+  4Trf,  TSi4TrfiT  -  4Trfi  T  Si  2TT(f2  “  fi  )T  -  4TTfi  T  Si  2iT(f2  +  f,  )T 
-  4TTf2  T  Si  2Tr(f2  +  fj  )T  +  4Trf2  T  Si  2TT(f2  —  fj  )t'1  (4.  22) 
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OzlT)  --T77^ 


^{iz  -  f 


f  I  1  1  1 

i)|^(fo+fi)  ^(fo+fz)  ^(fo  -  fi  )  -  t'o  ) 


,  cos  ZirTifo  -  f,  )  ,  cos  ZTTT(f:,  -  fn  )  .  cos  ZTTT(fo  +  fz  ) 
^(fo  -  fi  )  ^(fz  -  fo  )  ^(fo  +  fz) 

-  ^  +  iiTSiZTrT(fo  +  fz  )  -  nTSiZivTifo  +  fj  ) 


(4.Z3) 


+  ttT  Si  ZirT(fo  -  fj  )  +  wT  Si  ZirTifj 


-fo)| 


-  Z(f2  -f,  )  Z(fo+f,  )  Z(fo  +  f2  )  Z(fo  -  f,)  Z(f2  -  fo) 


cos 


ZTrT(fo  -  f,  )  .  cos  ZnT(f2  -  fo  )  ,  cos  ZTrT(fo  +  fz  ) 

iWo  -Xt^  ^  — mK-X)  '  ^ — w^nri 


(4.Z4) 


~~T(yJ ^  ^TTTCfo  +  fi  )  -  ttT  Si  ZTrT(fo  +  fj  ) 


+  TTTSiZTrT(fo  -£i  )  +  irTSiZirTlfz  -  fo  ) 

The  constant  y  is  equal  to  1.781.  The  function  Ci(x)  is  defined  as 

Ci(x)H  _  ^  £2^  du  (4.Z5) 

‘^x 

As  X  becomes  large  the  function  Ci(x)  converges  rapidly  to  zero  as 
(sinx)/x;  for  x  >  1  0,  Ci(x)  _<  0.1.  The  function  is  treated  fully  in 
Reference  19,  pp.  1-9. 

As  in  the  previous  system  the  rate  at  which  E/N  increases 
depends  upon  a  complex  function  of  the  chopping  frequency,  band¬ 
width  and  integration  time.  However,  some  qualitative  remarks  may 
be  made.  For  very  small  values  of  T  the  signal-to-noise  ratio  goes 
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to  zero  as  ;  for  large  values  of  T,  it  increases  as  the  square  root 
of  T.  Between  these  limiting  values  no  general  characteristics  can 
be  obtained  since  the  function  behaves  complexly. 

Let  us  consider  the  effect  of  the  chopping  frequency  and  the 
power  spectrums  of  the  cell  noise  upon  (4.20).  The  three  components 
of  the  noise  in  the  output  of  the  integrator  stage  are  given  in  Equa¬ 
tions  (3.6)  through  (3.9).  The  first  two  components  in  (3.6)  are  cross 
terms  involving  the  target  and  cell  noise  signals.  The  power  spectrums 
of  each  of  these  terms  are  nearly  identical  to  those  in  Figure  5  a,  b,  and 

c;  the  only  changes  are  in  the  substitution  of  e  i  or  e,  for  e  '  and  P 

'  “  *  ‘  1  max 

or  P  for  P  .  These  power  spectrums  are  shown  in  Figures 

2  max  max  r  t  o 

8  and  9.  The  third  component  of  (3.6)  has  a  power  spectrum  as 
shown  in  Figure  10.  The  power  spectrum  of  AN(t)  is,  of  course, 
the  sum  of  the  three  component  spectrums. 

From  the  above  three  figures  some  conclusions  can  be  reached 
for  optimum  correlator  efficiency:  1)  the  chopping  frequency  should  be 
kept  very  near  the  upper  band-pass  limit;  that  is,  fg  *  f2  :  and  2)  the 
chopping  frequency  should  be  as  high  as  possible.  These  conclusions 
follow  from  the  fact  that  the  integrator  stage  operates  as  a  low-pass 
filter  (see  (2.17)). 

If  the  signal-to-noise  ratios  ej  /vj  and  E2  /v^  are  assumed  to 
be  equal  (defined  as  c/v),  then  (4.20)  may  be  written  in  terms  of  the 
function  Q(T).  The  equation  becomes 


Q(T)  =  — 


sin  4'rTfo 
4TTfoT 


e 

V 


T 

sfoTT) 


(4.26) 


where  0(T)  is  given  by  (4.21)  through  4.24).  The  only  significant  dif¬ 
ference  between  this  equation  and  (4.12)  is  that  an  additional  factor 
e/v  has  appeared;  Q(T)  is  now  dependent  upon  the  original  signal-to- 
noise  ratio.  If  the  ratio  e  /  v  is  set  equal  to  unity,  then  a  graph  of  (4.26) 
may  be  made  similar  to  that  in  Figure  6.  The  graph  is  shown  in 
Figure  11.  The  signal  frequency  fg  has  again  been  given  the  same 
value  as  £2  .  As  in  Figure  6  the  function  Q(T)  increases  as  the  square 
root  of  the  time -bandwidth  product,  for  Q(T)  >  1.  However,  it  should 
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be  remembered  that  the  original  signal  -  to -nois  e  ratio  c  /  v  is  in¬ 
versely  proportional  to  the  square  root  of  the  bandwidth. 

A  graph  of  Q*  (T)  for  the  two-detector  c ros s -cor relation  system 
is  shown  in  Figure  IZ.  The  signal-to-noise  ratio  t/v  has  been  set 
equal  to  unity;  and  fg  ,  to  £2-  As  in  synchronous  rectification  the  improve¬ 
ment  in  signal-to-noise  ratio  for  Q(T)  >  1  is  nearly  independent  of  where 
the  bandwidth  is  located  in  the  frequency  spectrum.  That  is,  no  signifi¬ 
cant  advantage  is  gained  by  going  to  high  chopping  frequencies. 

Several  conclusions  can  be  obtained  by  comparing  Figures  6  and 
11.  The  curves  in  Figure  11  show  a  constant  displacement  downward 
of  a  factor  of  2.45  from  the  curves  in  Figure  6.  The  cause  of  these 
lower  Q(T)  values  are  the  additional  two  terms  in  the  noise  output  of 
the  correlator  (see  (3.6)).  A  second  conclusion  relates  to  the  dependence 
of  Q(T)  in  two-detector  cross-correlation  upon  the  original  signal-to- 
noise  ratio.  If  another  value  for  e/v  were  selected  for  Figure  11th 
curves  would  be  displaced.  The  trends  in  the  displacement  will  be  the 
following:  1 )  for  e  /  v  «  1 ,  Q  (T)  is  proportional  to  the  first  power  of  e  /  v; 
and  2)  for  tlv»  1,Q(T)  is  independent  of  the  original  signal-to-noise 
ratio.  Thus  synchronous  rectification  and  two-detector  cross -correla¬ 
tion  are  nearly  equivalent  for  e/v  »  1.  For  e/v  «  1»  the  improvement 
in  signal-to-noise  ratio  in  cross -correlation  is  less  than  that  in  syn¬ 
chronous  rectification  by  a  factor  of  2/(e/v). 
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V.  RETICLE  AND  SCANNING  SYSTEMS 

An  important  topic  to  be  considered  is  the  application  of  corre¬ 
lation  techniques  to  reticle  and  scanning  systems  performing  search 
and  tracking  operations.  The  primary  function  of  these  systems  is  to 
indicate  the  presence  and  position  of  the  target  within  a  specified  field 
of  view.  This  section  discusses  the  transfer  of  positional  information 
by  the  synchronous  rectification  and  two-detectqr  cross -correlation 
systems . 

In  reticle  systems  the  position  of  the  target  is  obtained  from 
information  inherent  in  the  signal  waveform.  The  sources  of  informa¬ 
tion  are  the  phase,  frequency,  and  amplitude  of  the  signal.  Various 
combinations  of  these  sources  are  used  to  locate  the  target  in  the 
field  of  view.  For  instance,  in  the  Sidewinder  system  the  phase  and 
amplitude  of  the  signal  provide  the  angular  and  radial  information 
necessary  for  locating  the  target.  When  applying  correlation  techniques 
to  reticle  systems,  the  immediate  problem  is  how  to  transfer  positional 
information  in  the  target  signal  and  yet  retain  a  relatively  high  informa¬ 
tion  rate.  Since  positional  information  depends  upon  the  shape  of  the 
signal  waveform  it  necessarily  follows  that  the  correlation  time  must 
be  short  compared  with  the  period  of  the  signal.  Otherwise,  the  origi¬ 
nal  shape  of  the  target  signal  is  lost.  The  difficulty  in  obtaining  signifi¬ 
cant  improvements  in  signal-to-noise  ratio  for  short  correlation  times 
(  -  10'^  sec)  can  be  easily  seen  from  Figures  6  and  11.  Only  for  sys¬ 
tems  having  wide  noise  bandwidths  (>  10  kc)  is  this  possible;  for  nar¬ 
rower  bandwidths,  the  improvements  are  very  small.  Improvements 
in  signal-to-noise  ratio  are  accompanied  by  a  reduction  in  information 
rate.  We  may,  therefore,  make  the  following  conclusion:  in  reticle 
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systems  requiring  high  rates  of  position  information,  it  would  appear 
that  correlation  techniques  will  find  no  application,  unless  the  condition 
(Af/fo)  ^  10^  is  met. 

The  scanning  systems  referred  to  are  those  using  the  scan¬ 
ning  spot  and  multiple -scanning  spot.  The  single  spot  obtains  target 
position  information  by  systematically  scanning  and  sampling  small 
elements  within  the  total  field  of  view.  The  small  element  repre¬ 
senting  the  instantaneous  field  of  view  localizes  the  target  in  the 
scanning  pattern.  The  multiple- scanning  spot  is  a  linear  array  of 
single  spots  performing  a  one -dimensional  scan  to  cover  the  field 
of  view.  Several  detectors  are  used  in  this  system,  each  one  cover¬ 
ing  a  portion  of  the  scanning  pattern.  In  both  of  these  methods,  the 
required  improvement  in  signal-to-noise  ratio  and  noise  bandwidth 
of  the  system  will  determine  the  minimum  dwell  time  of  the  scan¬ 
ning  operation.  The  dwell  time  is  to  be  equated  with  the  correlation 
time.  The  design  of  the  scanning  system  will  then  consist  of  a  compro¬ 
mise  between  the  parameters  of  information  rate,  total  field  of  view, 
and  instantaneous  field  of  view.  If  low  rates  of  target  position  informa¬ 
tion  are  acceptable  in  the  design  of  a  scanning  system,  then  application 
of  correlation  techniques  should  give  significant  improvement  in  the 
signal-to-noise  ratio  of  the  system. 

The  multiple  scanning-spot  method  implies  the  use  of  more 
than  one  correlator  in  the  scanning  operation.  This,  in  turn,  raises 
the  question  of  the  size  and  weight  of  correlation  systems.  Early  in 
i960  work  was  begun  on  the  production  of  a  miniaturized  correlator. 

The  result  was  a  correlator  composed  of  a  Hall-effect  multiplier  and 
an  RC  filter  acting  as  the  integrator.  A  photograph  of  the  system  is 
shown  in  Figure  13.  More  than  just  the  correlator  is  shown  in  the 
photograph.  The  multiplier  unit  is  located  in  the  upper  right-hand 
corner;  the  integrator  unit  is  in  the  lower  right-hand  corner.  The 
elements  on  the  left  side  of  the  photograph  are  transitor  amplifiers 
which  drive  the  Hall-effect  multiplier.  A  free-running  multivibrator 
circuit  coupled  with  a  relay  was  used  to  control  the  integration  time 
over  a  period  of  about  0.Z5  seconds  to  5  seconds.  The  scale  at  the 
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bottom  of  the  picture  is  in  inches;  the  third  dimension  of  the  correlator 
is  about  two  inches.  Figure  14  shows  a  circuit  diagram  of  the  correlator 
components  and  transistor  amplifiers.  The  frequency  response  of  the  mul¬ 
tiplier  is  limited  only  by  the  characteristics  of  the  coil,  since  the  res¬ 
ponse  of  the  Hall  element  itself  extends  well  into  the  megacycle  region. 

In  summary  we  may  conclude  that  correlators  may  be  obtained  which 
occupy  a  volume  of  four  or  five  cubic  inches  and  have  a  weight  of  a  few 
ounces . 

One  additional  application  of  correlation  techniques  should  be 
mentioned.  A  correlation  system  may  be  used  to  replace  the  cooling 
system  of  certain  detectors;  that  is,  detectors  whose  signal-producing 
mechanism  is  not  impaired  at  room  temperature.  However,  it  should 
be  emphasized  that  the  use  of  a  correlator  will  considerably  increase 
the  response  time  of  the  detection  system.  Recently,  a  successful  ex¬ 
periment  of  this  type  was  carried  out  with  photoconductive  InSb  cells 
operating  at  room  temperature.  With  a  noise  bandwidth  of  one  kilocycle 
and  an  initial  signal-to-noise  ratio  of  unity,  an  order  of  magnitude  im¬ 
provement  in  cell  detectivity  may  be  obtained  in  a  correlation  time  of 
0.2  0  seconds.  An  order  of  magnitude  improvement  in  detectivity  is 
equivalent  to  cooling  the  cell  to  dry-ice  temperature  (195®K).  See 
Reference  20,  p.  1467  f. 
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VI.  THE  EXPERIMENTAL  STUDY 

Many  correlators  have  been  constructed  for  the  purpose  of  uti¬ 
lizing  the  properties  of  the  correlation  function.  Among  these  may  be 
distinguished  two  classes:  the  sampling  correlator  and  the  continuous 
correlator.  The  first  class  of  correlators  operates  by  sampling  each 
time  dependent  function  entering  the  multiplier  stage  over  a  given 
period  of  time;  then  the  corresponding  samples  ^re  multiplied  and  added 

together.  Correlators  of  this  type  include  those  of  Cowley  (1949), 

22  23 

Singleton  (1950),  and  Cheatham  (1951),  The  second  class  of  corre¬ 
lators  performs  a  continuous  operation  of  multiplying  the  two  functions 

24 

together  and  then  integrating  (averaging)  the  result.  Stevens  (1950), 
Fano^^(1950),  Faran  and  Hills^^(1952),  and  John8on^^(1952)  have 
described  correlators  of  this  type.  The  correlator  developed  at  NSL 
belongs  to  the  second  class. 

The  two-detector  cross -correlation  system,  incorporating  the 
NSL  correlator  is  shown  in  schematic  form  in  Figure  15  and  in  pic¬ 
torial  form  in  Figure  16.  The  radiation  source  is  a  Sylvania  glow 
modulator  tube.  Type  R1131C,  modulated  by  a  sinusoidal  waveform. 

The  modulated  radiation  source  is  a  convenient  component  for  the  sys¬ 
tem  since  it  eliminates  the  necessity  of  a  chopper  and  provides  a  target 
signal  of  known  form.  The  two  PbS  cells  are  mounted  adjacently  at  the 
end  of  a  brass  tube.  Since  ac  pickup  is  quite  large  in  a  high-gain  network, 
the  radiation  source  and  detector  assembly  components  are  placed  within 
a  metal  container  (Figure  17).  The  preamplifier  units  following  the  PbS 
cells  are  Tektronix,  Type  122,  low-noise  amplifiers.  Bandwidth 
settings  include  the  followirg:  low  frequency  cutoffs  (3  db),  0.2,  0.8,  8, 
and  80  cps;  high  frequency  cutoffs  (3  db),  50,  250,  1000,  10000,  and 
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40000  cps.  Gain  settings  are  approximately  100  and  1000.  The  multi¬ 
plier  stage  of  the  correlator  is  a  Philbrick,  Model  HM,  combined  mul¬ 
tiplier/operational  fivefold.  This  hybrid  component  is  designed  for 
special  purpose  applications  in  analogue  computing  operations,  and 
consists  of  a  multiplier  unit  and  five  K2-W  operational  amplifier  units. 
The  multiplier  has  a  dynamic  range  of  —  50  to  +  50  volts  dc,  a  maxi¬ 
mum  d-c  error  below  one  per  cent  of  full  scale,  and  a  band  pass  of 
dc  to  15  kilocycles.  The  integrator  stage  of  the  correlator  uses  one 
of  the  operational  amplifier  units,  the  remaining  four  being  used  as 
additional  amplification  stages  beyond  that  of  the  preamplifiers. 

Oscilloscope  pictures  have  been  taken  which  illustrate  in  a 
graphical  manner  the  operation  of  the  correlator.  Figure  18a  shows 
a  picture  of  the  waveforms  entering  the  two  input  channels  of  the  cor¬ 
relator.  The  waveforms  consist  of  a  noise  component  alone:  no  radi¬ 
ation  (target)  signal  is  present.  Figure  18b  shows  the  output  of  the 
correlator,  the  integrator  unit  integrating  successively  over  a  period 
of  0,4  sec.;  that  is,  the  integrator  is  shorted  out  every  0.4  sec.  The 
square  divisions  along  the  abscissa  are  calibrated  to  O.Z  sec  per  division. 
The  zero  line  is  located  just  below  the  central  line  of  the  picture.  Figure 
19a  and  b  shows  the  change  in  the  output  of  the  correlator  when  a  small 
sinusoidal  target  signal  is  introduced  into  the  two  input  channels.  Des¬ 
pite  the  small  signal-to-noise  ratio,  the  output  of  the  correlator  (Figure 
19b)  shows  in  a  remarkable  way  the  presence  of  the  signal  component. 

The  successive  integrations  are  again  carried  out  over  a  period  of  about 
0.4  sec.  The  scale  along  the  abscissa  is  0.2  sec  per  division,  and  the 
zero  line  is  located  three  square  divisions  from  the  bottom  of  the  pic¬ 
ture. 

To  confirm  the  theoretical  analysis  of  the  two-detector  cross - 
correlation  system,  (3.19)  must  be  evaluated.  Since  the  target  signal 
is  sinusoidal,  the  only  functions  needed  to  evaluate  (3.19)  are  the  ampli¬ 
tude-normalized  power  spectrums  of  the  noise  entering  the  correlator. 
These  are  obtained  by  means  of  a  Br’uel  and  Kjaer  AF  spectrometer 
(see  Figure  20a,  b).  Since  Pi  (f)  and  P2  (f)  are  not  analytic,  a  computer 
program  was  set  up  to  evaluate  (3.12),  (3.14),  and  (3.17),  The  program 
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was  carried  out  on  an  IBM  704  computer  at  the  M.  LT.  Computation 
Center.  The  results  are  shown  in  Figures  21  and  22  as  the  solid 
lines.  The  signal  frequency  fo  was  500  cps,  the  initial  signal-to- 
noise  ratios  were  1.0  and  0.50,  and  the  equivalent  bandwidths  of  the 
noise,  Afj*  and  Af2*  ,  were  l.l6kc  and  1.04  kc.  The  signal-to-noise 
ratios  entering  both  channels  of  the  correlator  were  equal. 

From  the  experimental  study,  data  points  were  obtained 
giving  the  improvement  in  the  signal-to-noise  ratio  Q(T)  in  terms 
of  the  noise  bandwidths,  integration  time,  and  initial  signal-to-noise 
ratio.  Two  such  series  of  data  points  are  given  in  Figures  21  and  22. 

The  low  initial  signal-to-noise  ratios  of  1.00  and  0.50  were 
produced  by  noting  the  linearity  of  the  glow  modulator,  detector,  and 
preamplifier  system.  The  voltage  across  the  glow  modulator  tube  is 
proportional  to  the  preamplifier  output  voltage  over  a  wide  range  of 
voltages.  In  addition  to  the  initial  signal-to-no^se  ratio  e/v,  the  two 
quantities  required  to  determine  Q(T)  are  E(T)  and  N(T).  A  statistical 
sampling  procedure  was  followed  in  measuring  these  quantities,  the 
samples  being  recorded  on  a  Honeywell  Visicorder.  From  the  tapes, 
average  and  root-mean-square  values  of  E(T)  and  N(T),  respectively, 
were  obtained. 

We  proceed  now  to  an  error  analysis  of  the  experimental  data 
points  determining  the  function  Q(T).  Since  Q{T)  is  defined  by  the 
equation 


it  follows  that 


(6.1) 


AQ(T)  _  AE  .  Av  AN  Ae 

^  E  V  ~T 


The  root-mean-square  fractional  error  in  Q(T)  is  therefore 
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The  terms  on  the  right-hand  side  of  (6.3)  can  be  obtained  by  esti¬ 
mating  the  accuracy  of  the  measurement  of  each  term.  The  frac¬ 
tional  error  in  measuring  e  and  v  was  about  O.OZ  and  0.05,  re¬ 
spectively.  The  number  of  samples  taken  in  measuring  E(T)  and 
N(T)  was  Z5.  In  the  case  of  E(T)  the  error  varies  as  a  function  of 
the  correlation  time,  and  has  a  range  of  values  from  about  0.05  to 
about  0.03.  The  fractional  error  in  N(T)  is  approximately  equal 
to  the  reciprocal  of  the  square  root  of  the  number  of  samples  taken. 
With  a  sampling  number  of  Z5,  the  rms  error  is  about  0.20.  Thus 
the  principle  contribution  to  the  error  in  Q(T)  is  due  to  the  function 
N(T).  Using  the  above  estimates  (6.3)  can  be  evaluated  to  give 

A  K  21  per  cent  (6. 

If  the  scatter  of  the  data  points  around  the  computed  value  of 
Q(T)  is  assumed  to  follow  a  Gaussian  probability  distribution, 
then  68  per  cent  of  the  data  points  (or  2  out  of  3)  should  lie  within 
one  standard  deviation.  Figures  21  and  22  agree  with  this  assump¬ 
tion. 

We  may  conclude,  therefore,  that  within  the  limits  of  ex¬ 
perimental  error,  substantial  agreement  was  obtained  between 
the  theoretical  analysis  and  the  experimental  study. 
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VII.  SUMMARY 

In  this  report  we  have  reviewed  the  possibility  of  applying 
correlation  techniques  to  improve  the  signal-to-noise  ratio  of  in¬ 
frared  detection  systems.  General  expressions  have  been  derived 
which  describe  the  operation  of  two  correlation  systems  that  can 
be  used  for  this  purpose.  A  detailed  analysis  for  systems  having 
band-limited  white  noise  and  sinusoidal  target  signals  has  also 
been  presented.  The  most  promising  applications  of  these  tech¬ 
niques  appear  to  be  with  scanning  systems  and  mosaics  where 
high  information  rates  are  not  required.  An  experimental  study 
of  the  two-detector  cross -correlation  system  has  confirmed  the 
theoretical  analysis  of  this  system. 
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Dependence  of  Q{T)  upon  the  Position  of  the  Bandwidth  in  the  Frequency 
Spectrum  (Synchronous  Rectification) 
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Figure  8.  Power  Spectrum  of  ANj  (t) 
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Figure  13.  A  Miniaturized  Correlator  Using  a 

Hall-Effect  Multiplier  and  RC  Filter 
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igure  16.  The  Experimental  System  Shown  in  Pictorial  Form 
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Figure  17.  Radiation  Source  and  Detector  Assembly 
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(a)  Noise  Components  Entering  Both  Channels  of  the  Correlator 


(b)  Correlated  Noise  for  Integration  Time  of  about  0.4  sec 

Figure  18 
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(a)  Noise  Waveforms  with  Small  Signal  Component  Ent 


(b)  Correlator  Output  Showing  Presence  of  Signal. 


Figure  19 
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Figure  21  .  Experimental  Data  Points  and  Theoretical  Curve  for 
(€/v)  =  1.0 
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Figure  22.  Experimental  Data  Points  and  Theoretical  Curve  for 
(€/v)  =0.5 
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